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AN ANALYSIS OF THE EFFECTS OF AEROELASTICITY ON STATIC LONGITUDINAL STABILITY 

AND CONTROL OF A SWEPT-WING AIRPLANE 1 

By Richard B. Skooo 


SUMMARY 


A theoretical analysis has been made of the, effects of aero- 
elasticity on the stick-fixed static longitudinal stability and 
elevator angle required for balance of an airplane. The analysis 
is based on the familiar stability equation expressing the con- 
tribution of wing and tail to longitudinal stability. Effects of 
wing, tail, and fuselage flexibility are considered. Calculated 
effects are shown for a swept-wing bomber . of relatively high 
flexibility. 

Although large changes in stability due to certain parameters 
are indicated for the example airplane, the over-all stability 
change was quite small, compared to the individual e feels, due 
to the counterbalancing of wing and tail contributions. The 
efect of flexibility on longitudinal control for the example 
airplane was found to be of little real importance in Ig flight, 
although in turning flight the efect was found to be commensu- 
rate with the stability loss. 

INTRODUCTION 

In the past, airplane configurations and operating speeds 
have been such that prediction of longitudinal-stability and 
-control characteristics usually could be handled without 
regard to aeroelastic. effects. With higher flight speeds and 
the use of swept wings, aeroelastic effects are sufficiently 
important to pose major problems in airplane design. 
Although flexibility of swept wings has introduced the major 
problems to date, the flexibility of fuselage and swept-tail 
surfaces may introduce problems approaching equal impor- 
tance as speeds continue to increase. Much of the initial 
published work on predicting these effects appears to have 
been done by the British (refs. 1, 2, 3, 4, and 5). 

The object of the present study was to determine the 
magnitude of the stability loss likely to be encountered on 
swept-wing aircraft "of conventional configuration. To gain 
this end, it was decided to evaluate the aeroelastic effects 
for an airplane of relatively high flexibility. It was decided, 
further, to employ a simple method of stability analysis in 
order to seek a physical appreciation of the factors involved 
in the net stability change for the airplane. 

The results of the aforementioned study are presented in 
this report together with the method of analysis employed. 
The net stability change is shown together with the individual 
contributions due to flexibility of wing, tail, and fuselage. 


both including and neglecting the effect of inertial loads. 
The method of analysis is based on the familiar' stability 
equation expressing the contribution of wing and tail to 
longitudinal stability. The reader interested solely in the 
calculated results can turn directly to the section titled 
“Application to Example Airplane.” 
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NOTATION 

. b 2 

wmg aspect ratio? ^ 

ratio of net aerodynamic force along the airplane 
Z axis (positive when directed upward) to the 
weight of the airplane 

aerodynamic-center position, percent c behind 
leading edge 

wing span measured normal to the plane of 
symmetry, ft 

section chord parallel to the plane of symmetry, ft 

average section chord parallel to the plane of 
symmetry, ft 

rw 2 

J c?dy 

mean aerodynamic chord of wing, > 

Jo Cdy 

section lift coefficient 

lift coefficient, 

’ qS 

tail contribution to lift coefficient, 

qtS t 

wing lift-curve slope, per radian 

tail lift-curve slope, per radian 

pitching-moment coefficient, positive for nose-up 

, moment 
moment, 

pitching-moment coefficient at zero lift 

rate of change of pitching-moment coefficient with 
lift coefficient 


’ ft 


change in static stability parameter 


dC L 


due to 


flexibility (positive when stability is decreased 
by flexibility) 


• Supersedes recently desdassliled NAOA ESI A61C19 entitled “An Analysts of the Effectr of Aeroelastldty on Static Longitudinal Stability and Control of a Swept-Back-Wlng 
Airplane,” by Richard B. Stoog, 1951. 
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s, 

V 


tail incidence, i.,(X lr 4 z )+constant, (positive in 
same sense as a and also' relative to thrust axis 
of the rigid airplane), radians 
rate of change of tail incidence with wing angle of 
attack at constant A z (due to fuselage bending 
under the-aerodynamic load, L t ) 
rate of change of tail incidence with wing angle of 
attack at constant L t (due to fuselage bending 
under the loads imposed by the reaction of 
rear portion of fuselage and empennage masses 
to normal acceleration) 

rate of change of tail incidence with tail angle of 
attack at constant A z (due to fuselage bending 
under the aerodynamic load, L,) 
fuselage structural influence coefficient expressing 
change in tail incidence per unit normal ac- 
celeration, radians/# 

fuselage structural influence coefficient expressing 
change in tail incidence per unit tail load, 
radians/lb 

aerodynamic load on horizontal tail (positive when 
directed upward), lb . 

tail length (from airplane center of gravity to 
aerodynamic center of horizontal tail), ft 
free-stream dynamic pressure, lb/sq ft 
tail dynamic pressure, lb/sq ft 
wing area (including portion covered by fuselage), 
sq ft 

horizontal-tail area, sq ft 
true airspeed, ft/sec 


V tail volume, 

’ cS 

17 airplane weight, lb 

x distance from wing aerodynamic center to airplane 

center of gravity (positive when measured for- 
ward of center of gravity), ft 

y spanwise coordinate porpendicular to plane of 

symmetry, ft _ 

a angle of attack of wing root chord, radians 

a t angle of attack of tail root chord, radians 





elevator-effectiveness parameter, 


da, 

df. 


lift effectiveness of elevator, per deg 


dC m 

d8. 

8 < 

6 

de 

da 

6 


pitching effectiveness of elevator, per deg 

elevator angle (positive for down deflection),' deg 
downwash angle at the tail, radians 

rate of change of downwash angle at the tail with 
angle of attack 

pitching velocity, radians/sec 


V 

Aij 


nondimensional spanwise coordinate, 
spanwise shift in aerodynamic-center position for 
each wing panel, ^ 


Ac / 4 sweep angle of wing quarter-chord line (positive 

for sweepback), deg 

SUBSCRIPTS 

R rigid airplane 

F flexible airplane 

FUNDAMENTAL CONSIDERATIONS 

The 'purpose of this section of the report is to present tho' 
method of analysis used in obtaining the calculated results for 
the example airplane. The material is presented in three 
main subsections. In the first subsection, titled “Stability 
Equation,” the familiar stability equation expressing the 
contribution of wing and tail to static longitudinal stability 
is presented in a modified form to include factors which 
account for the airplane flexibility. In the second sub- 
section, titled “Effect of Flexibility on the Stability Param- 
eters,” methods are indicated and references are given to 
aid in evaluating the effects of flexibility on tho factors 
appearing in the stability equation. In 'the third subsebtion 
the effect of flexibility on longitudinal control is discussed. 

STABILITY EQUATION 

As is usually done, tho present report considers the index 
of static longitudinal stability to be the partial derivative 
dC„/dC L . In free flight, changes in lift coefficient at con- 
stant forward speed must always involve a curved flight 
path and hence must always involve pitcliing velocity. 
Further, in tho case of a flexible airplane, the normal acceler- 
ation associated with a curved flight path -will introduce 
deformations due to the loads imposed by the reaction of 
point masses to normal acceleration. It should be noted 
that, in general, the effects of such deformations will be in 
opposition to the effects associated with defoimations caused 
by increase in the aerodynamic loads only. In particular, 
the over-all effect of the mass reactions referred to is to pro- 
duce loads acting normal to the airplane plan' form which fie 
distributed over the wing and tail spans and along tho fuse- 
lage. In the discussion to follow, for brevity, these mass 
effects are referred to as inertial effects but should not bo 
confused with the effect of body inertia on dynamic longitu- 
dinal stability where the inertia of tho airplano as a whole is 
considered. The sign convention for the present analysis is 
shown in figure 1 . 



Figure 1 . — Schematic diagram defining positive directions' for the 
notation of the present analysis. 
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For ft rigid airplane, the stability equation expressing the 
contribution of wing and tail to bO n fdC L is usually written as 

The additional effect of pitching velocity on O m results in 

a (S$Hk®0. t ?s£ (2 > 

If equations (1) and (2) are modified to include effects of 
wing, tail, and fuselage flexibility and are then combined, 
the following equation can be written (see Appendix A) 

(S)J44 («o, KSs) J 4 ®t 

(3) 


As can be seen from these equations, flexibility of the air- 
plane structure is assumed to affect the following parameters : 

•C 

1 . Wing aerodynamic center position = 

c 

2. Wing lift-curve slope C La 

3. Tail lift-curve slope C Lai 

de 

4. Rate of chango of downwash at the tail ^ 

Also, additional parameters are introduced because of fuse- 
lage flexibility. Of the parameters listed, items 1, 2, and 4 
are influenced by wing flexibility with item 3 being influenced 
by tail flexibility. The effect of flexibility on these param- 
eters is discussed in detail later in this report. 

It is of some interest to note that the third term of the 
right-hand side of equation (3) and also the additional 
parameter (di,/da) L[ ‘m the second term of that equation do 
not appear in the evaluation of stability for a model in a 
wind tunnel where the model is physically restrained by the 
support system from developing any accelerations or angular 
velocities. However, since the former term is usually of 
small magnitude in comparison with the others in equation 
(3), it will be neglected in the present analysis. It should 
be noted, also, that the inertia of wing and tail has not been 
neglected in developing equation (3), although only the 
effect of fuselage inertia is shown explicitly. Actually, the 
direct effects of wing and tail inertia are considered' to be 
contained implicitly in (Ox, a ) F and (Xz, a J F - Details of this 

modification also will be referred to later in this report. 


aerodynamic center will cause a decrease in longitudinal 
stability. The geometric relation between any given span- 
wise shift and the associated chordwise shift along the mean 
aerodynamic chord is given by the following relation : 

4 (f)~ . (4 > 

It is apparent from this expression that the aerodynamic- 
center shift along the chord is most severe at high aspect 
ratios and high sweep angles for a given spanwise shift. 

A method for determining the span load distribution of a 
flexible sweptback wing, from which the aerodynamic- 
center position of the flexible wing can be determined, is 
given in reference 6. The shift in aerodynamic-center 
position for the tail is neglected since it only affects the 
value of l t used in the calculation of V. 

Wing lift-curve slope. — On a sweptback wing, the lift-curve 
slope for the flexible wing is usually less than that for a rigid 
wing. As a result, the angle of attack required to reach a 
given lift coefficient is higher for a flexible wing than for a 
rigid wing. For this reason, the tail also is subjected to a 
higher angle of attack (neglecting downwash considerations) 
so that the stability contribution of a rigid tail will be higher 
on an airplane with a flexible wing than on an airplane with 
a rigid wing. Consequently, the effect of reduction of wing 
lift-curve slope on the stability contribution of the tail is 
such as to increase the stability for an airplane with swept- " 
back wings. 

A method for determining the lift-curve slope for a flexible 
sweptback wing also is given in reference 6. 

Bate of change of downwash at the tail (de/da). — The rate 
of change of downwash at the tail is dependent upon the span 
load distribution associated with changing angle of attack 
so that changes in span load distribution due to flexibility 
may have an influence on the average downwash at the tail. 
The rate of change of downwash at the tail also depends 
upon the lift-curve slope of the wing, however, so that the 
over-all effect of wing flexibility on be/da cannot be stated 
even qualitatively for a general case without further analysis. 

Methods for predicting the downwash in the plane of the 
vortex sheet for low lift coefficients for any arbitrary continu- 
ous span load distribution are given in references 7 and 8. 
Although the value of befda may vary considerably in a 
direction perpendicular to the vortex sheet, it was felt that 
application of such methods to an estimation of changes 
in downwash was at least approximately correct. 

In the method of reference 7 (which is believed to be the 
simpler), the downwash variation of befdCz is calculated 
from a relation which is essentially 


EFFECT OF FLEXIBILITY ON THE STABILITY PARAMETERS 

Wing aerodynamic center. — On a sweptback wing of even 
moderate sweep, wing bending exerts the predominant 
influence on the aerodynamic span load distribution, resulting 
in an inboard shift in center of load for each wing panel. 
Duo to this inboard shift in center of load, the aerodynamic- 
center position for a flexible wing will be ahead of that for 
the rigid wing. As is well known, any forward shift in 
480101—08 10 


de 1 / CiC \ 

ZC L 2 (5) 

where the terms a-,* are influence coefficients given in reference 
7, and ( [c t c/C L Cat) n are loading coefficients corresponding to 
given spanwise stations as obtained from the span load 
distribution for the flexible wing. The span load distribu- 
tion for the flexible wing can be found from the method of 
reference 6. 
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To determine the stability parameter de/da, equation (5) 
is modified by introducing (C L ^j p . The resulting expression 
is 


de _ ( g Of ^ ( c,c \ 

bcT 2 A \C L c a h 


( 6 ) 


Tail lift-curve slope. — The effect of flexibihty on the 
lift-curve slope of a swept tail shows the same qualitative 
effects as for a swept wing. On a sweptback tail, the effect 
of aeroelasticity will he to cause a reduction in lift-curve 
slope so that the tail lift (and, consequently, the pitching 
moment) will be reduced at a given angle of attack. The 
effect of reduction in tail lift-curve slope on the stability 
of the airplane, therefore, is to reduce the stability, whereas 
the effect of reduction in wing lift-curve slope was shown 
to increase the stability. 

Tail lift-curve slope can be determined by the same method 
as for the wing as given in reference 6. 

Change in tail incidence due to fuselage bending. — From 
equation (3) it can he seen that fuselage flexibility affects 
the stability in general by the introduction of three additional 
terms giving the change in tail incidence due to fuselage 
bending. These terms are ( bi,[ba) Lt , and 

(bi,fda t ) Az . Expressions for these terms are developed in 
Appendix A but are shown below for convenience. 

/diA [1 - (£>e/da)„] ( 0 Lq( ) S t (Zi ,/di,)^ 

\teh~ (7) 

/SiA ^A 2 ) Lt [(C La ) F q/(W/S)] 

VWx,, l-(C Lat )^ t S t (bi t /bL t ) A2 W 

/cHA ( c z at ) F q,S,( ? >i,/dL ‘ )A * m 

\teJ Az l-(C Lat ) F q t S t (di t /?)L t ) Az W 


EFFECT OF FLEXIBILITY ON LONGITUDINAL CONTROL 

In the preceding sections of this report, the change in static 
longitudinal stability due to aeroelasticity has been discussed 
in some detail. It is also of interest, however, to consider 
the effect of flexibility on elevator angle required for balance 
(A m =0). There are two effects which must be considered: 

1. The primary effect of the stability changes just dis- 
cussed. 

2. The secondary effect of elevator deflection in the 
absence of any stability changes. 

The second effect is introduced by the deflection of hori- 
zontal tail and fuselage under the load produced by elevator 
deflection. 

The first effect can be evaluated graphically by plotting 
the stability curve (C m vs. C L ) for the rigid airplane and also 
the family of stability curves for the flexible airplane as 
calculated at several values of dyn ami c, pressure. In turning 
flight, the elevator angle for balance, then, can be calculated 
from the individual stability curves directly. In straight 
flight, where the dynamic pressure changes with lift coeffi- 
cient, the elevator angle variation can be found from the 
variation of C„ with C L given by plotting the flight variation 
of C L with q across the family of curves for the flexible 
airplane. 


The second effect can be evaluated by means of the follow- 
ing equation derived in Appendix B. 


5 e =6, 


(bC L JbS,) R r _/ djA 
(bC Lt /t>s.) F L l \mJa 


(PO# 8 '_ 


( 10 ) 


In this equation 5 a represents the elevator angle required 
(with flexibility of the horizontal tail and fuselago present) 
to maintain the same lift on the tail as that produced by an 
arbitrary elevator deflection (o, 0 ) on the rigid airplane. The 
ratio of elevator effectiveness for the rigid tail (pC L J<58,) n 
to the elevator effectiveness for the flexible tail ,) F 

can be found by the relation which is derived in Appendix C. 


(pC L JbS e ) F , AG;,, qt 

pc L jt ) «.)* c L ' R i+H 


where k and A C^JCl^ are constants -with magnitude de- 
pendent on the structural rigidity of the tail. Evaluation of 
k and AC L ' t jCz. tR (defined in Appendix C) involves a knowl- 
edge of the span load distribution due to elevator deflection 
and the span load distribution due to symmetric twist dis- 
tribution. Span load distributions for wings and flaps of 
arbitrary plan form are given in reference 7 for symmetric 
flap deflection. Span load distributions due to symmetric, 
continuous twist distribution are given in reference 9. The 
calculation procedure is similar to that contained in refer- 
ence 6. 

The procedure used to calculate the elevator angle required 
for balance may be summarized as follows : 

(1) Determine the stability curve for the rigid airplane for 
a given center-of-gravity position. 

(2) Using the stability given by (1) and the stability 
change given by equation (3), determine the slope of the sta- 
bility curve for the flexible airplane. Determine the change 
in (7,^ by integrating the aeroelastic loadings caused by 
built-in twist and camber as found by the method of reference 
6. Note that the effect of dead weight on C„ 0 is accounted 
for by using stability curves which include the effects of 
inertia. With the slope and intercept thus determined, plot 
stability curves for the flexible airplane at several values of 
dynamic pressure. 

(3) Plot the flight variation of C L with q across the family 
of curves so obtained. 

(4) Deter min e the elevator anglo required for balance 5 (() 
corresponding to the variation in C m with C L so determined. 

(5) Obtain the final elevator angle required for balance by 
means of equation (10). 


APPLICATION TO EXAMPLE AIRPLANE 

The method of analysis indicated in the preceding section 
of the report has been applied to an example swept-wing 
airplane known to have a relatively flexible structure. Com- 
pressibility considerations in regard to the effect on span 
load distribution were neglected in evaluating the aeroelastic 
effects since a preliminary estimate showed them to bo of 
second order (compared to the primary effect of dynamio 
pressure) for the particular configuration studied. The 
ratio of tail dynamic pressure to free-stream dynamio pres- 
sure was assumed to be equal to 1.0. The airplane con- 
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Figure 2. — Geometric characteristics of example airplane. 


figuration is shown in figure 2 with the pertinent geometric 
parameters indicated. The sweep angles of wing and tail 
are 35° and 33°, respectively; the wing aspect ratio is 9.43, 
wing taper ratio 0.42; tail aspect ratio 4.06, tail taper ratio 
0.423; and tail volume 0.672. The effect of the engine 
nacelles on the aerodynamic span load distribution was 
neglected as was the effect of the fuselage. The elastic axis 
for the wing is located at the 38-percent-chord line and for the 
tail is located at the 50-percent-chord line. The variation 
of pertinent structural characteristics across the semispan of 
of wing and horizontal tail is presented in figure 3. The 
structural influence coefficients associated with fuselage 
flexibility are defined by the following data: 


(I?) =-°- 0000342 > deg/to 

(&Xr 0A5 ’ deg/9 


EFFECT OF FLEXIBILITY ON STABILITY 

In the discussion to follow, the material is presented in 
the following order: first, an evaluation of the effects of 
flexibility on each of the parameters appearing in the static- 
longitudinal-stability equation, and second, an evaluation of 
the over-all effect of flexibility for the airplane as a whole. 
The individual effects are presented in figures 4 through 9 
and are summarized in figure 10. The over-all effect is 
presented in figure 11. 

The effect of changes in the stability parameters on the 
stability of the example airplane has been evaluated by 
assuming that only the parameter under consideration is 
affected by flexibility. In this way, the effect of changes in 
each parameter can be assessed individually. In the discus- 
sion of each parameter which follows, effects due to the 
action of aerodynamic loads only are considered first, fol- 
lowed by consideration of the modifying influence of inertia. 


Wing aerodynamic center. — The shift in wing aerodynamic- 
center position due to wing flexibility is shown in figure 4 
together with the associated stability change. Inasmuch as 
the discussion is being restricted for the moment to effects 
due to aerodynamic loads only, merely the curves shown for 
the weightless wing need be considered, since these curves 
are for zero inertial effect. As can be seen, the aerodynamic 
center with inertia absent moves forward from the rigid- 
wing value of 32 percent of the mean aerod 3 Tiamic chord 
until at a dyn ami c pressure of about 800 pounds per square 
foot the aerodynamic center is almost at the leading edge 
of the mean aerodynamic chord. As can be seen from the 
values of stability change, the effect of aerodynamic-center 
shift in itself is very large. For example, at a dynamic 
pressure of 500 pounds per square foot, the neutral point of 
the wing has shifted forward by 20-percent chord, which of 
itself would introduce a serious stability problem. 




0 .2 .4 .6 .8 1.0 

Spanwise station, y 

(a) Wing. ' 


Figure 3. — Variation in pertinent structural characteristics across the 
semispan for the wing and horizontal tail of the example airplane. 






130 


REPORT 12g8 — NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 




JD 



AT 

„g» 

> 

“O 

s 

Q 



(b) Horizontal tail. 
Figure 3. — Concluded. 


It is of some interest to know how much of the stability 
change with inertia absent is due to bending deformations 
and how much is due to torsidnal deformations. In order to 
aid in this comparison, a curve has been included in figure 4 
showing the contribution of torsion alone. As can be seen, 
torsional deflections are stabilizing, while bending deflections 
are destabilizing. The contribution due to torsion is seen 
to be much smaller than that due to bending. The relative 
importance of torsion and bending, however, depends on the 
ratio of torsional to bending rigidities and location of the 
elastic axis, and hence would not necessarily be the same for 
all airplanes. An equally important factor to consider is 
the effect of sweep angle. The extremes of zero sweep and 
,90° sweep best illustrate the point, since for zero sweep only 
torsion is a factor, while for 90° sweep only bending is a 
factor. 

The effect of inertia on the location of the aerodynamic 
center for the example wing also is included in figure 4 for 
airplane wing loadings of 70 pounds per square foot and 100 
pounds per square foot. As can be seen from the figure, the 
effect of wing inertia is only mildly alleviating. Although 
the relieving effect in this case is shown to be rather small, 



Dynamic pressure, q, Ib/sq ft 


Figure 4. — Variation of wing aerodynamio center with dynamic pres- 
sure for the example airplane. 


the effect for other airplanes may not be of similar magnitude, 
since the inertia effect depends upon the ratio of wing weight 
to total airplane weight in addition to the Bpanwiso weight 
distribution previously mentioned. By reference once more 
to the case of a tailless airplane, it would appeal that wing 
inertia would have a much greater relieving effect in that 
case, since more of the total airplane weight is in the wings 
than for conventional airplanes. 

Wing -lift-curve slope. — The ratio of flexible to rigid wing 
lift-curve- slope and the associated increase in tail contribu- 
tion to longitudinal stability is presented in figui'6 5 as a 
function. - of dynamic pressure';’ As Before; the curves for 
the weightless wing represent the caie ofztero inertia effect. 
At a- dynamic pressure of 500 p (hinds pel Square! ‘foot the 
lift-curve slope with inertia absent is reduced to 64 percent 
of the rigid-wing value. The associated increase in • tail 
stability contribution amounts to a rearward neutral-point 
shift of 25 percent. At this same dynamic pressure, the 
stability contribution of the wing aerodynamic-conter shift 
(with inertia absent) was shown to be a forward shift of 20 
percent, or almost the same magnitude, so that the two whig 
factors so far discussed would appear to bo largely canceling. 
Whether canceling of these effects will exist in general for 
all configurations cannot be determined at this time. Calcu- 
lations for a fighter configuration of markedly different 
geometric and structural characteristics, however, resulted 
in essentially the same relation between these wing factors. 
An interesting extreme to consider is the case of tho tailless 
airplane for which the second term of the static-longitudinal- 
stability equation does not exist. In this case no canceling 
of these effects will be possible so that the net stability cliango 
will be due solely to wing aerodynamic-center shift. 

Since the effect on stability of reduction in wing lift-curvo 
slope is large with inertia absent, it is of interest to consider 
the relative contribution of bending and torsion. For this 
purpose a curve is included in figure 5 showing tho contri- 
bution of torsion alone. As can bo seen, the contribution 
of torsion to the lift-curve slope is an increase, while the 
larger effect of bending is a decrease. Tho associated 





ANALYSIS OF THE EFFECTS OF AEROELASTICITY ON STABILITY AND CONTROL OF A, SWEPT-WING AIRPLANE 131 


stability changes are shown to be a decrease due to torsion 
and a much larger increase due to bending. 

The effect of inertia on the wing lift-curve slope also is 
included in figure 5 for airplane wing loadings of 70 and 100 
pounds per square foot. As can be seen, the effect of wing 
inertia is much the same as was stated for aerodynamic 
center (i. e., mildly alleviating). 

Rate of change of downwash at the tail. — The variation 
along the swept-tail span of the rate of change of downwash 
in the plane of the vortex sheet with inertial effects absent is 
presented in figure 6 for several values of dynamic pressure. 
As indicated earlier in this report, the curves are based on a 
method which is applicable only at low lift coefficients. The 
location of the tip of the horizontal tail is indicated in the 
figure. As can be seen, large changes in downwash are 
indicated behind the outer sections of the wing and in the 
plane of symmetry; however, the average downwash over 
the tail is changed only slightly. The change in average 
downwash depends on the ratio of tail span to wing span. 
The downwash factor, 1— (del da), based on the average 
downwash over the tail is presented in figure 7 as a function 
of dynamic pressure along with the associated change in 
stability contribution of the tail. As can be seen from the 
figure, the change in downwash factor is very slight, being 
of the order of 5 percent at the highest dynamic pressure 
considered. The stability change, as would be expected, is 
correspondingly small and relatively unimportant compared 
to the other stability factors so far discussed. 

The effect of inertia on downwash at the tail is not shown, 
since the relatively larger downwash changes associated 
with the aeroolastic effects due to aerodynamic loads only 
were shown to be unimportant. 

Tail lift-curve slope. — The ratio of flexible to rigid tail 
lift-curve slope and the associated decrease in tail stability 
contribution is shown in figure 8 with similar curves for the 
wing shown for comparison. As can be seen from the figure, 
the effect of flexibility on tail lift-curve slope is not so pro- 
nounced as that for the wing, and, as a consequence, the 
effect on the stability contribution of the tail is also corre- 
spondingly loss. As can also be seen from the figure, the 
inertial effect on the tail is similar to that for the wing. 8 

Tail incidence change due to fuselage bending. — The 
variation with dynamic pressure of an over-all stability term 
expressing the change in tail incidence due to fuselage bend- 
ing is shown in figure 9. The curves showing the contribu- 
tion of aerodynamic loads only are indicated by (di,/ba) i( ==0. 
Curves are presented considering the effect of fuselage 
flexibility alone and also including the effects of wing and 
tail flexibility. For comparison, curves of average down- 
wash and stability change due to downwash change are also 
presented. At the higher values of dynamic pressure the 
fuselage factor becomes of the same order of magnitude as 

* Tho figure also Indicates the extent to which the Inertial effect varies with wing loading, 
since the curve for the weightless wing could Jost a3 wen be labeled IF/S— to . x physical 
explanation here Is that on airplane of Infinite mass will experience rero At under lifting loads, 
so that tho effect of flexibility for 1175= » 13 thg same as If the airplane were physically re- 
strained and merely pivoted about the center of gravity. An airplane having finite 
however, will experience a normal acceleration so that the wing Inertial loads will affect the 
wing deflection and hence Cl ^ 
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Figubb 5. — Variation of the ratio of flexible to rigid wing lift-curve 
slope with dynamic pressure for the example airplane. 



Fionas 6. — Variation along the swept-tail span of the rate of change cf 
downwash in the plane of the vortex sheet for several values cf 
dynamic pressure for the example airplane as affected by aero- 
dynamic loads only. 



Figube 7. — Variation of the downwash factor [1 — (de/da) a d and the 
associated stability change with dynamic pressure for the example 
airplane as affected by aerodynamic loads only. 
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Dynamic pressure, q, Ib/sq ft 


Figure 8. — Variation of the ratio of flexible to rigid tail lift-curve 
slope and associated stability change with dynamic pressure for the 
example airplane, with similar curves for the wing shown for com- 
parison. 


the average rate of change of downwash at the tail and there- 
fore is seen to be of comparable importance. The effect of 
including wing and tail flexibility in the fuselage factor is 
to lower the factor slightly as shown. By reference to the 
stability curves, it can be seen that the stability change due 
to fuselage bending is of much greater importance than that 
due to downwash change, as would be expected from the 
comparison shown on the upper part of the figure. It can 
also be seen that the effect of wing and tail flexibility is to 
alleviate the stability decrease due to fuselage flexibility. 

As can be seen from the figure, the effect of inertia on the 
fuselage factor is very large and consequently of considerable 
importance. It will be remembered that the effect of inertia 
on the wing and tail factors was only slight by comparison. 
It is interesting to note that consideration of inertia and of 
all the flexibilities involved results (for the example airplane) 
in a fuselage factor equal essentially to zero even though the 
aerodynamic contribution is large. In these estimates of 
inertial effects, the influence on fuselage factor of wing and 
tail inertia has been neglected, since these effects are of higher 
order for this airplane. 

Summary of effects. —The effects of wing, tail, and fuselage 
flexibility on the longitudinal-stability index i>C„fdC L of the 
example airplane are summarized in figure 10, showing the 
important individual effects which have been discussed. The 
upper set of curves presents the effects due to aerodynamic 
load only while the lower set of curves also includes the effects 
due to inertial loads for an airplane wing loading of 70 pounds 
per square foot. As can be seen from the figure, all the effects 
are destabilizing except the effect of reduction in wing lift- 
curve slope on the stability contribution of the tail. The 
stability changes due to wing aerodynamic-center shift and 
reduction in wing lift-curve slope are shown to be by far the 
largest effects of those shown. Both these results are shown 
to be true whether inertial effects are included or not. 

The over-all effect of flexibility on the static-stability index 
bC m fdC L for the example airplane is shown in figure 11. 
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Figure 9. — Variation with dynamic pressure of the factors associated 
with fuselage flexibility and the stability change duo to theso factors 
for the example airplane. 



Dynomic pressure, q, Ib/sq ft 


Fioube 10. — Summary of the individual effects of tho various param- 
eters involved on stability of the example airplane as affected by 
aerodynamio loads only and also as affeoted by both aerodynamic 
and inertial loads. 


The curves presented show the changes due to aerodynamic 
loads only and also include the effects of inertial loads for 
an airplane wing loading of 70 pounds per square foot. Duo 
to the nature of the second term of the stability equation, 
the effects shown in figure lOUre not all additive algebraically; 
therefore, the curves of figure 11 were obtained by allowing 
all the factors in the equation to vary simultaneously. As 
can be seen from the figure, inertial effects reduce con- 
siderably the stability change which would otherwise occur 
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Figure 11. — Over-all effect of flexibility on the static-stability index 
dC m /bCL for the example airplane as affected by aerodynamic loads 
only and also by both aerodynamic and inertial loads. 

for this airplane. For example, the stability change at a 
dynamic pressure of 500 pounds per square foot is reduced, 
due to inertia, from a neutral-point shift of about 17 percent 
to about 7 percent. When the large individual effects are 
remembered, the over-all change due to flexibility for the 
airplane (as shown by the inertial curve) is seen to be rela- 
tively small. 

It is of interest to consider the effect of flexibility on the 
flight test index b8 e fbC h (which is associated with bC m fbC L ), 
sinco 3 

/S5A (bC m fbC L ) F 
\bCrJr (bCJbS.) F 

Consideration of the above relation alone leads to the pos- 
sibility of counter-balancing between the effects of aero- 
elasticity on bC m fbC L and on bC m /b5 e -, therefore, figure 12 is 
presented to show the over-all change in bS,fbC L as a function 
of q for an assumed static margin for the rigid airplane of 
0.10. Curves are presented for bC„lb8 e = (dO^/dS,)* and for 
bC m /b 8 e = (bC m /b8 e )f to show the influence of that parameter 
on A(b8 e /bC L ). As can be seen from the figure, the variation 
of A(b8,/bC L ) with q for the example airplane is governed 
primarily by the stability change rather than by loss in 
control effectiveness. That this conclusion will be true, in 
general, for airplanes which suffer loss in stability with 
increase in q perhaps is not immediately obvious, but can 
be seen quite readily from recognition of the fact that the 
intersection of the two curves shown in figure 12 at a dynamic 
pressure of about 600 pounds per square foot is due to the 
flexible airplane having become neutrally stable. It is 
apparent that the parameter bC„fb8, can have no effect 
on bSgfbC L when an airplane has neutral stability since 
bd'/bC^O. 

EFFECT OF FLEXIBILITY ON LONGITUDINAL CONTROL 

The effect of flexibility on the elevator angle required for 
balance in lg flight has been evaluated for the example air- 
plane for a static margin of the rigid airplane of 0.10. The 
semigrapliical analysis described earlier was used. The case 

> Tho parameter (oCmfdCOr Is given by equation (3) of this report. The parameter 
(0 C m fb}.)r con be found from equation (10) since, from the nature of the analysis leading to 
that equation, It Is dear that 

St 


of turning flight is not considered here since the results for 
that case are believed to be sufficiently summarized in figure 

12 just discussed. The variation of elevator angle with 
dynamic pressure for steady level flight is presented in figure 

13 (a) for the case where up-elevator only is required and 
in figure 13 (b) for the case where both up- and down- 
elevator are required. 4 Curves are shown for both the rigid 
and flexible airplane. As can be seen from the figure, the 
effect of aeroelasticity for the case of up-elevator only is 
quite small. This results from the counteracting effects of 



Dynamic pressure, q, Ib/sq fl 

Figure 12. — OVer-all effect of flexibility on the rate of change of ele- 
vator angle with lift coefficient in turning flight for the example 
airplane assuming, static margin for the rigid airplane of 0.10. 



Dynamic pressure, q Ib/sq ft 

(a) With only up-elevator angle required (C„ o =0). 

(b) With both up- and down-elevator angle required (C„ o = +0.03). 

Figure 13. — Variation in elevator angle for balance with dynamic 
pressure in steady level flight for the example airplane when con- 
sidered to be both rigid and flexible assuming static margin for the 
rigid airplane of 0.10. 

flexibility on the slope of the stability curve, the change in 
due to camber and twist, and the loss in elevator effective- 
ness. For the case where both up- and down-elevator are 
required, the more positive value of causes a more 
downward increment of elevator deflection for the flexible 
case than for the rigid because of the loss in elevator effective- 
ness due to flexibility. 

* The determining factor In each case Is the value of for the rigid airplane. In case 1 

Cmg was assumed equal to lero, while for ease 2, was assumed equal to +0.03. In cal- 
culating the change in C ^ due to flexibility for both cases, values ol 1° washout and c. 0 = 
—0.010 (section pitching moment due to camber) were used. 
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The effect of flexibility of horizontal tail and fuselage in 
' modifying the elevator angle required for balance is shown in 
figure 14 as calculated for the example airplane. The effect 
shown in this figure is included in the curves of figure 13. 
The elevator-angle ratio increases almost linearly reaching 
a value of about. 1.9 at a dynamic pressure of 800 pounds per 
square foot. 

In all of these calculations, torsion of the elevator has been 
neglected; that is, the surface is assumed to be infinitely 
rigid to hinge moments. 
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Figure 14. — Variation of the elevator angle ratio 3«/5. 0 with dynamic 
pressure for the example airplane. 


izations of the results presented herein. Although the over- 
all aeroelastic effect on stability for the example airplano 
was found to be small compared to the individual effects, 
it cannot be said that like calculations for any airplano will 
also yield a small effect. It can be said, however, that, for 
any practical swept-wing airplane with a tail, the stability 
change due to a shift in wing aerodynamic center will bo 
destabilizing, while the change in wing lift-curve slope will 
be stabilizing, so that a certain amount of counterbalancing 
between these major effects will always be present. As can 
be seen from the simple stability equation employed in this 
analysis, the degree of completeness of the counterbalancing 
depends directly on the size, plan form, and location of the 
tail as they affect the stability contributed by the tail. 
Therefore, in the design of airplanes for which wing flexi- 
bility would be expected to exert a large influence on sta- 
bility, it would appear that a minimum over-all aeroelastic 
effect may be accomplished more advantageously by design 
changes to the horizontal tail than by such changes to the 
wing. 

The calculations presented herein with regard to the effect 
of flexibility on longitudinal control show the effect to be of 
little real importance for the example airplane in lg flight, 
although in turning flight the effect is shown to bo commen- 
surate with the stability loss. 


CONCLUDING REMARKS 

The effect of flexibility on static longitudinal stability and 
control of an airplane is such as to preclude any real general- 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Mar. 19, 1961 


APPENDIX A 


EVALUATION OF PARAMETERS ASSOCIATED WITH FUSELAGE FLEXIBILITY 


In order to show how the parameters associated with 
fuselage flexibility enter the longitudinal-stability equation, 
equation (3) of the text will be derived here. If the pitching- 
moment contribution of wing and tail are considered in 
nondimensional form (i. e., in terms of C L , C Lt , etc.), the 
pitching-moment coefficient of the combination can be 
written as: 

(CJr=«i) r (§)-(C e ,),Tf (Al) 

Since a flexible airplane is being considered, (0 L ) F and 
{C T ^ P should be expressed in terms of parameters applying 
to a flexible structure, so that 

(A2) 
(A3) 

(A4) 


where 


(pL tt ) r a 


a t — a F — er-Mt+Zt y- 


(Cj F =(C L y(=)- (0 Lat ) F (a F - eF +i t+ !£) (A5) 

Differentiating with respect to 0 L and factoring (G La ) F from 
the second term, equation (A5) becomes 


/dO,\ _/*\ ( c Q y r. /a a ^ &n 

wJrvA (GO, L 



Before proceeding further, it is necessary to evaluate 
To do this, the tail incidence (i,) can be written 


Az Li+ Qa z ) Li Az (A7) 


where (di t ji) Lt ) Az and (bi,fbA z ) Li are structural influence 
coefficients associated with fuselage bending. Differentiating 
equation (A7) with respect to a F , 


In equation (A4), it should be recognized that the downwash 
angle e F may be affected by the changes in span load distri- 
bution associated with flexibility and also that the tail 
incidence i, is affected by fuselage flexibility. Substituting 
equation (A4) into equation (A3) followed by substitution 
of the modified equation (A3) into equation (Al), 



It should be noted here that, although equation (A6) con- 
tains the partial derivative bi,fda F , the total derivative of 
equation (A7) must be taken because of the dependence of 
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e, i t , and 6 on a F in the expression for L t which can be written 
as 

L,= (C La ^ F (ctF — fjp+ii+p- 6^ gtS\ (A9) 

The derivative of equation (A9) -with .respect to a F is then 

,+k+vr] (A io > 

For notational consistency with equation (A7), cU/da F and 
dd/da F can be written as befba F and bdfba F , since at con- 
stant q the quantities e and 6 are functions of a F only. 
With these substitutions in equation (AlO) and knowing 
that 

dA z _{C La ) F q 

da F W/S 

from differentiation of 


(All) 


a (CM 
Az ~~W[S~ W/S 


(A12) 


equation (A8) can be written after some rearrangement as 


where 


di,_fbi t \ fbi,\ U bd . /bi t \ 

dcc F Vx*t/A z Vbcryba)^ 

/bi t \ [l-Wte) F ] (C Lat ) F q t S t &i,/bL t ) Az 
V b^) A]S ~ 1 - ( C Lai )q t S t (bi,/bL t )A, 


(A13) 


• (pL a ^ F <lt>J t \W,I UJ^ tUz 

\ (C Lat ) F q t S t (bi t /bL t ) Az 

)a x ~ 1 ~ (C^ F q t S l (bi,fbL t ) Az 

/diA (c X'/bA s ) Lt [(C La ) F q/(3V/S)] 

\ba) Lt 1 (bi,/bLi) Az 

By substitution of equation (A13) into equation (A6) and 
after rearrangement, a final expression for bC m fbC L may be 
written as 


bC n \ /x\ ( 0l «) f W fbi\ . fbi t \ ,fbit\ 

bC L )p \cJ F (C La ) p 1_ \d«/F + \ b«A z + \ ba) Lt \ q 


(A14) 


APPENDIX B 

EVALUATION OP ELEVATOR ANGLE REQUIRED FOR BALANCE 


In evaluating the elevator angle required for balance 
(C m — 0) for a flexible airplane, it is necessary to consider the 
effect that elevator deflection will have on aeroelastic dis- 
tortion of the horizontal tail and fuselage since such distor- 
tion will change the elevator angle required to balance the 
pitching moment existing with elevator neutral. If the 
elevator is assumed to be deflected by an arbitrary angle 5, 0 
with the fuselage and horizontal tail fixed in position, the 
lift coefficient on the tail due to elevator deflection may be 
written as 


C, 


-( b Jh) s 

bS, A ‘° 


(Bl) 


If the horizontal tail is now allowed to relax (with the fuselage 
still fixed in position), the lift on the tail will change due to 
distortion of the structure. The elevator angle required to 
maintain the same lift on the tail as that given by equation 
(Bl) is defined by 


so that 


Cw)>= a '-Cw)/ 

. » 

n * pc L jbsj F 


(B2) 


(B3) 


If the fuselage now is allowed to relax, the lift on the tail 
will change due to a change in tail incidence. The additional 
elevator angle required to maintain the lift at the value 
given by equation (Bl) can be written as 


A i t 


where 


13 K)r 


(B4> 


The final elevator angle required is then 

5«=5 4l +5, 2 


(B5) 


(B6) 


By substitution in equation (B6), the following expression is 
obtained for the elevator angle required on a flexible airplane 
to maintain the same lift on the tail produced by an arbi- 
trary elevator deflection on a rigid airplane 


* , r( dC V &5 «)s fbi,\ (bC L JbS e ) R q t Sr 

* ,0 L (bC Lt /bs.) F V>lJ Az k). 


Since 


( ) (bCiJ bS*) 

Ke) F (bC Lt [b a ‘) 


(B7) 


(B8) 


equation (B7) can be written more conveniently as 

<B9> 
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APPENDIX C 


EVALUATION OF LIFT EFFECTIVENESS FOR FLEXIBLE TAIL 


The lift effectiveness of an elevator (i>C L[ fdS t ) for the case 
of a flexible horizontal tail will differ from that for a rigid 
tail due to distortion of the structure . 1 The streamwise 
twist of the structure can be foimd in reference 6 by applying 
a relaxation approach to the problem. Using this approach, 
the twist distribution for a sweptback tail can be written 
in series form as 

e(v)=fo(v)~‘ Ae 1 (j 7 )+A£ 2 (ij)+ . . . (Cl) 

where. 

twist of the flexible wing 

twist produced by the loading for the rigid wing 
twist produced by the loading associated with ^(ij) 
twist produced by the loading associated with Aej(ij) 
etc. 

Since t(jj) can be written in series form, it is apparent that a 
similar expression can be written for the lift coefficient 
produced by a given elevator deflection for the flexible wing, 
so that 

Cz. t =C Lt -AC Lti +AC Lti . . . (C 2 ) 

where 

Cl. lift coefficient for the flexible tail 

<F 

Cl 1r lift coefficient for the rigid tail 

A Ct lift coefficient obtained by integrating the loading 

associated with < 3 ( 77 ) 

A Cl, lift coefficient obtained by integrating the loading 

associated with Aei(y) etc. 

If the terms of equation (Cl) aro related by a constant of 


1 Elevator distortion Is neglected In this analysis. 


«(*) 

<0(17) 

A«i(ij) 

AcaO?) 


proportionality (which is the usual case ), 2 equation (C 2 ) can 
be written as 


where 


CL'=Cr.,-AC Lli (l-K+K> . . .) 

A Cl, 


(03) 


K= 


A C L 


Equation (C3) simplifies to 


C L , =C L , - AC . i 


■‘i (i+a:) 


(04) 


since the series 1—K+K 2 ... is merely an expansion of 
1/1 +K. The ratio of the lift effectiveness for the flexiblo 
tail to that for the rigid tail can bo obtained from equation 
(C4) by dividing by C L[r . The equation becomes 


&CL t /?)5') p 0 Lif a c Lll 1 

&C L JbK) R C LtR AC LtR l+K 


( 05 ) 


Equation (C5) can be simplified since K and A Cl, are 
proportional to q c and becomes 

t /d5.) y A<7 S/ g « q, A<7 U, q, , 

,/&5.) fi Cl, b l+iK/qdq, Cl 1r 1 +kq, 

where A Cl\ and k are evaluated for unit dynamic pressure 
at the tail. 

1 The line of load application for the first term win) is the centroid of tho chorihvlse loading 
produced by elevator deflection while tho line of load application for the remaining terms is 
the aerodynamlo center of the section. Tho relative contribution of bending and torsion, 
therefore, is not the same in <o(v) ns in succeeding terms. Calculations made to dote, however, 
have shown tho constant of proportionality to apply to tho relation between ttM and A.iin) 
as well as to the terms beyond 4 o(e). There may be some configurations for which tho 
proportlonaUty will be limited to A, 1(7) and succeeding terms, however. 
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